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Abstract. - We show how the renormalized force correlator A(u), the function computed in the func- 
tional RG (FRG) field theory, can be measured directly in numerics and experiments on the dynamics of 
elastic manifolds in presence of pinning disorder. For equilibrium dynamics we recover the relation ob- 
tained recently in the statics between A(u) and a physical observable. Its extension to depinning reveals 
interesting relations to stick-slip models of avalanches used in dry friction and earthquake dynamics. The 
particle limit (d — 0) is solved for illustration: A (it) exhibits a cusp and differs from the statics. We 
propose that the FRG functions be measured in wetting and magnetic interfaces experiments. 



Models involving elastic objects driven through random media are important for numerous physi- 
cal systems and phenomena including magnets [1], superconductors [2], density waves [3], wetting [4], 
dry friction [5], dislocation and crack propagation [6], and earthquake dynamics [7]. There has been 
progress in qualitative understanding of, e.g. the existence of a depinning treshold for persistent motion 
at zero temperature T = 0, scale invariance at the threshold and the analogy with critical phenomena, 
collective pinning and roughness exponents, avalanche motion at T — 0, and ultra-slow thermally 
activated creep motion over diverging barriers. These phenomena are predicted by theory, i.e. phe- 
nomenological arguments [2], mean field models [8], functional renormalisation group [9-12], and 
were seen in numerical studies. Experimental evidence for creep motion was found in vortex lattices, 
in ferroelectrics, and in magnetic interfaces in ferromagnets [1,13]. Some cases exhibit clear discrep- 
ancies with the simplest theories, e.g. the depinning of the contact line of a fluid [4, 14]. Even when 
agreement exists, much remains to be done for a precise comparison. 

Recent theoretical progress makes these quantitative tests possible. For interfaces, powerful algo- 
rithms now allow to find the exact depinning threshold and critical configuration on a cylinder [15] 
and to study creep dynamics [16]. The functional RG has been extended beyond the lowest order (one 
loop), and it was shown that differences between statics and depinning become manifest only at two 
loops [11,12]. The FRG is the candidate for a field theoretic description of statics and depinning, be- 
yond mean field. It captures the complex glassy physics of numerous metastable states at the expense 
of introducing, rather than a single coupling as in standard critical phenomena, a function, A(it), of 
the displacement field u, which flows to a fixed point (FP) A*(u). This FP is non-analytic, as is the 
effective action of the theory. At a qualitative level, A(u) can be interpreted as the coarse-grained cor- 
relator of the random pinning force and its cusp singularity at the FP, A'(0 + ) = — A'(CP), is related to 
shock singularities in the coarse grained force landscape, reponsible for pinning. Until now however, 
comparison between experiments, numerics and FRG was mostly about critical exponents. 
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The aim of this paper is to make precise statements concerning the physics of dynamical FRG 
and propose experimental and numerical tests. Recently a relation was found [17] between the FRG 
coupling functions A(u) = —R"(u) and observables, suggesting a method to measure these functions 
in the statics. The idea is to add to the disorder a parabolic potential (i.e., a mass m) with a variable 
minimum location w. The resulting sample-dependent free energy V(w) defines a renormalized ran- 
dom potential whose second cumulant is proved to be the same R(w) function as defined in the replica 
field theory - deviations arising only in higher cumulants [17]. This holds for any internal dimension 
d of the elastic manifold, any number of components N of its displacement field u(x), and any T. At 
T = 0, the (minimum energy) configuration u(x; w) is unique and smoothly varying with w, except 
for a discrete set of shock positions where u(x; w) jumps between degenerate minima. The limit of a 
single particle in a random potential (d = 0) maps to decaying Burgers turbulence, and the statistics 
of the shocks can in some cases be obtained, yielding exact result [17] for A(u). 

This method was used recently [18] to compute numerically the zero-temperature FRG fixed-point 
function A(u) in the statics, for interfaces (TV = 1), using powerful exact minimization algorithms. 
Random bond, random field and periodic disorder were studied in various dimensions d = 0, 1, 2, 3. 
The results were found close to 1-loop predictions and deviations consistent with 2-loop FRG. A linear 
cusp was found in any d and the functional shocks leading to this cusp were seen. The cross-correlation 
for two copies of disorder was also obtained and compared to a recent FRG study of chaos [19]. The 
main assumptions and central results of the FRG for the statics were thus confirmed. It is important to 
extend these methods to the dynamics of pinned objects and to the depinning transition. 

In this Letter we extend the method of Ref. [17] to the dynamics. Using a slow, time-dependent, 
harmonic potential we show how the various terms in the effective dynamical action identify with 
the FRG functions. The T > equilibrium dynamics reduces to the same definition as used for 
the statics. We describe the extension to depinning at T = 0. There the manifold is pulled by a 
quasi-static harmonic force (i.e. a spring of strength noted m 2 ), and we show how the statistics of 
the resulting jumps directly yields the critical force and the FRG functions, and how they converge to 
fixed forms as m — ► 0. The model is similar to some stick-slip models used e.g. in dry friction [5,20] 
and earthquake dynamics [7]. The present method provides a different way to look at these problems 
in numerics and experiments, in addition to giving a precise meaning to quantities computed in the 
field theory. In particular we discuss the identification of the critical force, the statistics of the jumps, 
using for illustration a graphical construction in d = 0. There we compute the FP functions for each 
universality class. These exhibit a cusp which we find is rounded by a finite velocity. These effects 
could be tested in experiments, as discussed at the end. 

We consider the equation of motion for the overdamped dynamics of an elastic manifold parame- 
terized by its time-dependent displacement field u(x, t): 

r)d t u(x,t) = F x [u(t);w(t)} (1) 
F x [u; w] — m 2 (w — u(x)) + cV x u(x) + F(x, u(x)) 

where F x [u(t);w(t)] is the total force exerted on the manifold (we note u(t) = {u(x, t)} xe9 d the 
manifold configuration, x being its <i-dimensional internal coordinate); 77 is the friction coefficient and 
c the elastic constant. Here at the bare level, the random pinning force is F(x,u) = —d u V(x,u) 
and the random potential V has correlations V(0, x)V(u, x') — i?o(u)<5( d ) (x — x'). We consider first 
bare random bond disorder with a short-ranged Rq(u). At non-zero temperature one adds the thermal 
noise (£(x, t)£(x' , t')) = 2r]T5(t — t')5 d (x — x'). We have added a harmonic coupling to an external 
variable w(t), a given function of time (in most cases we choose it uniformly increasing in t). This is 
the simplest generalization of the statics, where w(t) = w is time-independant. It is useful to define 
the fixed-u; energy 

/2 
d d x^-{u{x)-wf + V(x,u(x)) (2) 
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associated to the force F x [u;w] — — ■ If w(t) is an increasing function of t the model represents 

an elastic manifold "pulled" by a spring. Quasi-static depinning is studied for dw/dt — > + . 

We first describe qualitatively how to measure the FRG functions and later justify why the relation 
is expected to be exact. Consider the observable w(t) — (u(t)), where u(t) = L~ d J d d x u(x, t) is the 
center of mass position, and (. . .} denotes thermal averages, i.e. the ground state at zero temperature. 
It is the shift between the translationally averaged displacement and the center of the well, i.e. the ex- 
tension of the spring. It is proportional to the pulling force on the manifold, hence to the translationally 
averaged pinning force minus the friction force, i.e. w{t) — u{t) = m~ 2 (r]v(t) — f F(x, u(x, t))) (if 
we use periodic boundary conditions inside the manifold). Of particular interest are: 

w(t) - (u(t)) = m- 2 f av (t) (3) 
Hi) - <«(*)>] Hi') - {u(t'))] c = m- 4 L- d D w (t,t') , 

where connected means w.r.t. the double average (...). If we consider a function w(t) such that 
dw(t)/dt > 0, one can also write: D w (t,t') — A w (w(t),w(t')). As written, the function A w may 
in general depend on the history w(t). However we expect that for fixed L, m and slow enough w(t), 
e.g. w(t) = vt with v — > + , one has A w (w(t),w(t'j) — > A(w(t) — w(t')). This function A(w — w'), 
which is independent of the process w(t), is the one defined in the F.T., as we will justify below. 

Let us start with non-zero temperature, T > 0, and consider a process w(t) so slow that the system 
(with a finite number of degrees of freedom (L/a) d ) remains in equilibrium. In practice it means that 
iMl -C u(L) where tz is the largest relaxation time of the system, and u(L) its width. The above 
definition is then consistent with the one from the statics, where it was shown that one can measure the 
equilibrium free energy in a harmonic well with fixed w (or its generalization to an arbitrary w(x)), 
defined through e~ v ( w ^ T = J T>[u] e _w »M/ T , and extract from it the pinning energy correlator 

R(w). This can be done by measuring the second cumulant [21] V(w)V(w') = R[w — w'], with 
R[w] = L d R(w) for a uniform parabola w(x) = w, and using that R — R [17]. One equivalently 
obtains the force correlator A(w) via the equilibrium fluctuations of the center of mass (u) w at fixed 
w, i.e. (w — (u) w )(w' — (u) w >) = m~ 4 L~ d A(w — w'). In the statics it is easy to show that A(w) 

/,' The potentiality of this function breaks down in the driven dynamics, or at depinning, as 
discussed below. 

Let us note at this stage that a second definition can be given using two "copies". Consider two evo- 
lutions u(x; wi) and u(x; W2) driven by two (slow) processes w\{t) = W2(t) +wof fixed separation, 
in the same disorder sample. Then define 

( Wl (t) - <th(t)»(t«2(i) - (u 2 (t))) C = m- 4 L- d A t (w) (4) 

which is now an equal-time correlation. For a slow equilibrated motion at T > 0, it identifies with the 
static definition. The general case is discussed below and in [24]. 

Let us now describe T — depinning, and restrict to TV = 1. Quasi-static depinning is studied as 
the limiting case where dw/dt — > + . One starts in a metastable state Uo(^) f° r a given w = wq, i.e. 
a zero-force state F x [uq (x);w) = which is a local minimum of H wo [u] with a positive barrier. One 
then increases w. For smooth short-scale disorder, the resulting deformation of u{x) is smooth. At 
w = w\, the barrier vanishes. For w = wf the manifold moves downward in energy until it is blocked 
again in a metastable state u\{x) which again is a local minimum of H Wl [it]. We are interested in the 
center of mass (i.e. translationnally averaged) displacement u = L~ d J d d x u(x). The above process 
defines a function u(w) which exhibits jumps at the set Wi. Note that time has disappeared: evolution 
is only used to find the next location. The first two cumulants 

w — u(w)=mr 2 f c , (w — u(w))(w' — u(w') — m~ A L~ d A{w — w') (5) 
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Fig. 1 - Main plot: Construction of u(w) in d — 0. The pinning force F(u) (bold black line). The two quasi-static 
motions driven to the right and to the left are indicated by red and green arrows, and exhibit jumps ("dynamical 
shocks"). The position of the shocks in the statics is shown, for comparison, from the Maxwell construction 
(equivalence of light blue and yellow areas, both bright in black and white). The critical force is 2/L times 
the area bounded by the hull of the construction. Inset: The same construction for the forward motion of the 
discretized model. 



allow a direct determination (and definition) of the averaged (m-dependent) critical force f c and of 
A(w), in analogy to the statics. Note that u(w) depends a priori on the initial condition and on its 
orbit but at fixed m one expects an averaging effect when w is moved over a large region. This is 
further discussed below. Note that the definition of the (finite-size) critical force is very delicate in the 
thermodynamic limit [22]. Here the quadratic well provides a clear way to obtain a stationary state. 

Elastic systems driven by a spring and stick-slip type motion were studied before, e.g. in the 
context of dry friction. The force fluctuations, and jump distribution were studied numerically for a 
string driven in a random potential [20]. However, the precise connection to quantities defined and 
computed in the field theory has to our knowledge not been made. The dependence on m for small m 
predicted by FRG, A(w) = m e ~ 2 '' A(wm~^) is consistent with observations of [20] but the resulting 
A (to) has not been measured. Fully connected mean-field models of depinning also reduce to a particle 
pulled by a spring, together with a self-consistency condition, around which one can expand [10]. As 
discussed below, our main remarks here are much more general, independent of any approximation 
scheme, and provide a rather simple and transparent way to attack the problem. 

For the qualitative discussion it is useful to study the model in d = 0, i.e. a particle with equation 
of motion 

rjdtu = m 2 {w — u) + F(u) . (6) 

In the quasi-static limit where w is increased slower than any other time-scale in the problem, the zero 
force condition F(u) = m 2 (u — w) determines u(w) for each w. The graphical construction of u(w) 
is well known from studies of dry friction [5]. When there are several roots one must follow the root 
as indicated in Fig. [2 where F(u) is plotted versus m 2 (u — w). This results in jumps and a different 
path for motion to the right and to the left. Let us call A the area of this hysteresis loop (the area of 
all colored/shaded regions in FigQ. It is the total work of the friction force when moving the center 
of the harmonic well quasi-statically once forth and back, i.e. the total dissipated energy. The above 
definition of the averaged critical force Q, assuming the landscape statistics to be translationally 
invariant (hence replacing disorder averages by translational ones over a large width M) gives 

27 ™ 2 f M , n A 

f c = TO (w-u w ) = — J dw (w - u w ) = — , (7) 

where we have used J u dw = J w du. One can check that for m — > this definition of f c becomes 
identical to the one on a cylinder, f d , which for a particle (d = 0) is 2f d = f£ — = max M F(u) — 
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min u F(u) with 2/dM = lim m ^o A(m). (Since A depends on the starting point, this definition 
holds after a second tour, where the maximum (minimal) pinning force was selected). Finally, one 
can compare with the definition of shocks in the statics. There, the effective potential is a continuous 
function of w. Therefore, when making a jump, the integral over the force must be zero, which 
amounts to the Maxwell-construction of figure^ 

One can compute f c and A(w) in d = for a discrete force landscape, Fj, independently dis- 
tributed with P(F), and i integer. u(w) is then integer and defined in the inset of figure^ The process 
admits a continuum limit for small m, which depends on the behaviour of P(F) in its tails (nega- 
tive tail for forward motion). One obtains [24] the distribution of u(w), P w (u)du = e~ au,( ' u ^da w (u) 

where a' w (u) = /_„ P(/)d/ and a w (— oo) = 0. One also obtains the joint distribution of 
(u(w),u(w')), P w;w ,(u,u') = (a'Ju) - fl ' K ,(H))«U«>-»-(»)-«»'l»')+«^(»)%' - u ) + S(u' - 
u)a' w ,(u)e~ au, ( u ) for w > w' . Define A(w) —: m 4 p^ l A(w/p m ) and f c =: f° + cm?p m . This yields 
two main classes of universal behaviour at small to. The first contains (i) exponential-like distributions 
with unbounded support i.e. lnP(/) wj_ > _ QO — A(— /) 7 (for which = ((\nm~ 2 )/A)~) and (ii) 
distributions with exponential behaviour near an edge P(f) ~ e- A ^ +fo ^e(f + f ) (with 7 < and 
fc = fo~ {(lnm- 2 )/A)5). For both (i) and (ii) the FP function is A (x) = ^ + Li 2 (l - e x ) + ^ 
and p rn = p7n := l/(|7|A 1 / 7 TO 2 (lnm _2 ) 1_ ~), c = 7^ the Euler constant. The first class has 
£ = 2 up to log-corrections. The second class contains power-law distributions near an edge P(f) = 
Aa(a-l)(f + f ) a - 2 6(f + f ),a> 1, for which c = -T(l + 1), m 2 p m = (to 2 /A)^ and / c = /„. 
The FP depends continuously on a with A (to) = ~r(l + i)T(l + A, w 7 ) + wT(l + A) e -™ T + 
J °° dye-to +, "^ + ^r(l + A, y7) w h ere T(a,x) = f x °° dzz^e-"; it has C = 2 - 2/a [23]. Hence, 
despite the fact that d = is dominated by extreme statistics (e.g. the distribution of p^(w — u(w)) 
converges to the Gumbel and Weibul distributions for class I and II respectively) it still exhibits some 
universality in cumulants and in all classes A (it) has a cusp non-analyticity at u = 0. We have checked 
the above scaling functions and amplitudes numerically, with excellent agreement. 

We now come back to the interface, d > 0, and note that the manifold in the harmonic well can be 
approximated by (L/L m ) d roughly independent pieces with L m ~ 1/m. The motion of each piece 
resembles the one of a particle, i.e. a d — model, but with a rescaled unit of distance in the u direction, 
u m ~ L^ n ~ to - ^. The "effective-force" landscape seen by each piece becomes uncorrelated on such 
distances, and its amplitude scales as F m ~ m 2 u m . Hence one is in a bulk regime not dominated 
by extremes, i.e. A(w) probes only motion over about one unit. It is easy to check on Fig. 1 that an 
arbitrary initial condition joins the common unique orbit after about one correlation length. Hence 
the d = model suggests that starting the quasi-static motion in uo and driving the manifold over 
w ~ L^ n should then result in all orbits converging. Hence the definitions and are equivalent 
for N = 1. An interesting crossover to d = behaviour and extremal statistics occurs if L < L m . 

Note that the averaged critical force, defined in 0, should, for d > 0, go to a finite limit, with 
f cirri) = f£° + Bm 2 ~'' from finite size scaling. Although f c is not universal and depends on short- 
scale details, one easily sees that —md m f c (m) depends only on one unknown scale. We note that 
the definition l|5} coincides with the one proposed recently as the maximum depinning force for all 
configurations having the same center of mass uq [22]. Since u — w is a fluctuating variable of order 
(L/L m ) _d / 2 , the definition is the same as the above in the limit where L — > 00, before m — > 0. The 
single w distribution is obtained from the distribution of w — u(w) if all modes have a mass. 

Measurements of A(w) reveal interesting features in any d. At the bare level, the disorder of 
the system is of random-bond type (i.e. potential). As the mass is decreased, one should observe a 
crossover from random-bond to random-field disorder. Also a finite velocity should round the cusp 
singularity. These features are well visible in the d = toy model as illustrated in Fig.|2k (quasi-static 
evolution in a model with F, ; = Vi+i — V% and Vi uncorrelated) and Fig.|2p (Langevin dynamics at 



6 



EUROPHYSICS LETTERS 




Fig. 2 - (a): Main plot. The measured A(u) for RB disorder distributed uniformly in [0, 1], and rescaled such that 
A(0) = 1 and dw|A(«)| = 1. From bottom (which has J °° duA(u) ~ 0) to top the mass decreases from 
m 2 — 1 to 0.5 to 0.003. One observes a clear crossover from RB to RF. 

(b): Inset. A' (if) for a driven interface at vanishing (blue, A'(0) < 0) and finite velocity (red, A'(0) = 0). 

finite v) and can be obtained analytically in that case [24]. 

We now sketch the exact relation to the FT definitions. From now on we use condensed notations 
u(x, t) = u xt and so on. The bare MSR action functional can be parameterized as S[u, u] = u ■ g~ x ■ 
u + u- [u] — i{t • f?( )[u] • u + 0(u 3 ) with g xy an arbitrary (time independent) symmetric matrix 
(standard choice being g~ 1 = q 2 + to 2 ), A^[u] xt = r\d t u xt and B^[u] xt , x >t' = 2riT8 XX '6tt> + 
Ao(u xt — u x i t ')5 xx >. We denote u-v := J xt u xt v x t (and additional index contraction for TV > 1), A^ 
and B^ are respectively vector and matrix functionals. The effective action T[u] can be parameterized 
identically with ^4(°)[it] — > A[u] and B^[u] — > B[u]. It is obtained from the generating function: 
W[w, w] = In J T>[u] T>[u] exp(— S[u, u] + u ■ g^ 1 ■ w + w ■ g^ 1 ■ u) through a Legendre transform: 
W[w, w] + T[u, u] = u ■ g^ 1 ■ w + ils ■ g^ 1 ■ u. It can be expanded [21] as: W[w, w] — w ■ g^ 1 ■ 
w — w ■ A[w] + i-w • B[w] ■ w + 0(w 3 ). This functional directly generates correlations (0 and 

(0, in a more general form: w xt - {u xt ) w = g xy A yt [w] and (w xt - (u xt ) w ){w x > t i - {u x < t ')w) C = 
9xygx'y'B y t,y't' [w]. Hence A and B are observables which can be measured, i.e. for a uniform w x t = 
m, fav(t) = jit J y A yt [w] and D w (t,t') = J yyl B yt .y>t>[w}, which for slow w t should go to f c 
and A(w t — w t ' ) respectively. The question is how to relate them to the effective action, i.e. generalize 
the relation R = R from the statics. To this aim we perform a Legendre transform. Details are 
given in [24]. The result is = A[w], (where we have defined u[w] :— w — g • A[w], i.e. 

Uxt[w] := (u x t)w) an d B[w] = (du/dwY ■ B[u[w]] ■ dii/dw with du/dw = 1 — (V^Afiy])' • g. Now 
consider w uniform in space w x t = w t . Then A yt [w] is y-independent. In the limit of infinitely slow 
monotonous wt one expects Ayt[u>] — > f c sgn(w) hence u x t[w] — > Wt + sgn(w)f c /m 2 (with usually 
f c = at T > 0). From STS it implies A[w] = A[w] in that limit for uniform w. Similarly one 
finds D w (t, t') — J2 xy B x t,yt' [w] for infinitely separated times (with fixed w t — w' t ) since du/dw — 

g Su g^ — > 1, i.e. the disorder-averaged response function becomes trivial at zero frequency (due again 
to STS). Hence the measured A(w) in (|3J is - to all orders - the one defined in the FT. 

We have generalized [24] the above method to a manifold driven in N dimensions (e.g. a flux line 
in a 3D superconductor). For a particle, and fixed m, we have seen numerically that different initial 
conditions converge. The Middleton theorem no longer holds, and particles can pass by each other. 
To probe transverse motion and correlations A(w) in the transverse direction, we use the two-copy 
definition ©. Finally thermal rounding of depinning, creep and crossover from statics to depinning 
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can be studied more precisely by this method. 

To conclude let us propose that A(w) be measured directly in experiments, which would represent 
an important test of the theory and the underlying assumptions. Creep and depinning of magnetic 
domains in thin films with surface step disorder have been investigated using imaging [1]: adding 
a magnetic field gradient should allow to confine the interface in an effective quadratic well, whose 
strength and position can be varied (hence probing both statics and dynamics). In contact lines of fluids 
it is capillarity and gravity which provide the quadratic well, and provided large scale inhomogeneities 
can be controlled, A(u) could be measured from statistics on lengths larger than the capillary length 
(i.e. L m here). We hope this will stimulate further numerical [25] and experimental [26] studies. 

We thank A. Fedorenko, E. Rolley, A. Rosso and S. Moulinet for useful discussions, the KITP for 
hospitality and support from ANR (05- BLAN-0099-01). 
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